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ANALOG COMPUTER SIMULATION OF A PARASITICALLY LOADED 
ROTATING ELECTRICAL POWER GENERATING SYSTEM 
by Leonard J. Gilbert 
Lewis Research Center 

SUMMARY 

An analog computer representation was developed to simulate an alternator with a 
parasitically loaded speed controller. The simulation permits the study of the dynamic 
and steady-state behavior of the average values of the electrical characteristics of the 
power generating system. A detailed simulation of the electrical performance of a para- 
sitically loaded power generator had not been developed previously. 

The discussion describes a balanced three-phase electrical generating system with 
balanced loads. The computer representations for all of the electrical components are 
presented. The identity of the components is retained in the simulation. Speed regula- 
tion, voltage regulation, load power factor, saturation effects, and amortisseur damp- 
ing are factors included in the representation. A good correspondence between computer 
and theoretical performance results, indicating that the simplifying assumptions are ac- 
ceptable. 

The simulation is designed specifically to accommodate the study of paralleled para- 
sitically loaded alternators. A transformation to direct and quadrature reference axes 
is used in the analysis. 


INTRODUCTION 

An analog computer simulation was developed to simulate an electrical generating 
system consisting of an alternator with a parasitically loading speed controller. The 
simulation described provides a means of studying the dynamic performance of such an 
electrical power system. The computer output signals are proportional to average val- 
ues of the electrical quantities simulated. Analog computer simulations of a parasiti- 
cally loaded alternator have been reported by Tew (ref. 1) and Holt (ref. 2). But the de 
tailed simulation described here has not been developed previously. In particular the 



simulation is designed to be useful for studying the performance of paralleled parasiti- 
cally loaded generating systems. 

Several power generating systems are being developed for space applications (refs. 
3 to 5). The systems operate with constant total electrical loads provided by means of 
parasitic loads. Parasitic loading is an operational mode resorted to when the energy 
in the system energy source is relatively unlimited. As a consequence of its unlimited 
energy, the source can be used to control the generating system as well as to supply the 
electrical requirements of a useful load. 

In these parasitically loaded equipments the parasitic load is used to complement 
the operational or useful load of the alternator so that the sum of all loads on the alter- 
nator is constant, as illustrated in figure 1. By thus limiting the load torque on the tur- 



Figure i - Schematic diagram of ideal load distribution of parasitically 
loaded alternator. 


bine drive, the turbine speed is also limited. The magnitude of the parasitic load re- 
quired to balance a particular useful load is determined by a controller which senses the 
electrical frequency of the alternator output. This frequency varies inversely with the 
useful load applied to the alternator. By appropriate design of the sensing and loading 
circuits, the parasitic load and its control is made to serve as a speed control. 

The parasitic loading control approach is resorted to in closed-cycle space power 
generating systems in order to avoid rapid control of source energy and precise flow 
control of the working fluid in the system. Such controls are necessary if the working 
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fluid flow is used to regulate the turbine speed. However, high operating temperatures, 
the fast reaction speeds required, and the corrosive nature of some working fluids make 
it undesirable to vary flows continuously for turbine regulation. 

In a series of studies for the U. S. Air Force in 1957 the Energy Conversion Group 
at the Massachusetts Institute of Technology developed a simulation of aircraft electrical 
power generating systems (ref. 6). The approach used here is modeled after that devel- 
opment. However, in the aircraft systems simulations parasitic loading of the alternator 
is not a factor, whereas it is an integral part of some space power systems. 

The discussion which follows describes a balanced three-phase alternator with bal- 
anced loads. The computer representations for all of the electrical subsystem compo- 
nents are presented, followed by a sample set of displays of characteristic electrical 
signals from the simulation. 

The specific system simulated is comprised of the following principal components: 

(1) Turbine driven alternator 

(2) Voltage regulator and exciter 

(3) Parasitic load speed controller 

(4) Useful load 

The identification of each component of the system described is retained in the simula- 
tion. Physical parameters which determine the characteristics of the system compo- 
nents are identifiable. Speed regulation, voltage regulation, load power factor, satura- 
tion effects, amortisseur damping, and saliency are factors included in the representa- 
tion. The simulation is flexible enough so that it can be augmented to include the 
dynamics of the drive, and startup and shutdown operation, although this is not done 
here. 


MATHEMATICAL REPRESENTATION OF ELECTRICAL SYSTEM COMPONENTS 

System Configuration 

The program described simulates the electrical system illustrated in figure 2. A 
three-phase synchronous alternator is turbine driven with constant power. The electri- 
cal output of the alternator is applied to (1) a load representing the desired output, 
termed the useful load, and to (2) a parallel load, termed the parasitic load. Because 
of the user's varying demands, the useful load varies. The parasitic load is used to 
keep the total load constant by balancing it with the input. By keeping the total load 
power constant the parasitic load functions as a speed control. A voltage regulator- 
exciter energizes the alternator field and regulates the terminal voltage of the alterna- 
tor as the volt-ampere demand varies. 
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Figure 2. - Electrical power generating system block diagram. 


Analytical Approach 

In order to study the dynamic performance of the electrical quantities of the system 
the set of simultaneous differential equations which describe the system must be solved. 
The solution of these equations is extremely complicated unless use is made of a trans- 
formation suggested by Blondel's "two-reaction method" and developed by Park (ref. 7). 
The equations described in this report have been simplified by resolving the variables 
into components along direct and quadrature axes. This is a standard analytical ap- 
proach, and discussions of this approach are widely available in modern literature on 
synchronous machines (ref. 5). However, the use of quadrature reference axes for the 
variables is a significant feature of this simulation, and the feature which differentiates 
it from previous simulations of parasitically loaded generating systems. 

The simulation was undertaken for use ultimately in studying the dynamic behavior 
of paralleled parasitically loaded alternators. For paralleled systems information with 
respect to the relation between voltage and current and rotor is necessary for a com- 
plete understanding of performance. The use of the dq transformation not only greatly 
simplifies the solution of the system equations, but also facilitiates the determination of 
voltage- current phase relations. 

When used for the study of two paralleled systems, duplicate equipments described 
here are required, in addition to transformation and load sharing networks. 


Assumptions 

The following assumptions have been made to facilitate and simplify the analysis: 
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(1) There is sinusoidal air gap flux distribution and sinusoidal winding distribution; 
that is, the magnetomotive forces in the armature are sinusoidally distributed. 

(2) Phase loads are balanced. 

(3) Eddy current and hysteresis losses can be neglected. 

(4) Saturation effects take place only in the direct axis circuits. 

(5) Rotor damping can be represented by including two windings, one in the direct 
axis and one in the quadrature axis. 

(6) Friction and windage losses are neglected. 

(7) Turbine output power is constant. 

Assumptions (1) to (3) are basic assumptions made by R. H. Park (ref. 7) in the devel- 
opment of his synchronous machine analysis. That analysis provides the basic circuit 
equations used to describe the alternator. Assumption (4) reduces the complexity of the 
simulation. This reduction is desirable and, in view of the alternator performance data 
usually available, is warranted. The usefulness of the simulation is not substantially re- 
duced by this simplification. Assumption (5) is made because it is not possible to define 
a current path through a solid rotor. A fictitious set of damper windings is assumed. 

The windings are represented in the analysis by one winding in each of the quadrature 
axes. 


Alternator 


The conventional generator equations, in terms of direct and quadrature components, 
as developed by R. H. Park, are listed in this section. They describe the functioning of 
the circuits shown in figure 3. These equations are expressed in a per unit system 
chosen to make all the mutual reactances in any axis, d or q equal (ref. 8). (All sym- 
bols are defined in the appendix. ) 


v f = R f i f + pi// f 


v d = - R akd + p ^d - 
v q = - Vtq + P ^q + ^ d 
0 = fWkd + p ^kd 

0 = R kq\q + P ^kq 
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Quadrature axis 



Figure 3. - Schematic diagram of synchronous generator with two damper wind- 
ings. 


The flux linkages in the preceding voltage equations are as follows: 

^f = Vd + ^fL = Vd + X fL*f 

~ Vd + ^dL -l ^md ” X alAd 

^q = Vq + ^qL = ^mq ” X aL*tq 
^kd ~ ^md + ^kdL = ^md + X kd*kd 
^kq ~ ^mq + ^kqL ~ ^mq + X kq*kq 


The mutual flux linkages are 


Vd X md^f *td + *kd^ 
Vq = X mq^ _i tq + *kq^ 
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The active and reactive power generated at the terminals of the alternator and the 
electromagnetic torque developed are given by 

P t = Vtd + Vtq 

Qfc = "Vtq + v q l td 

T l = ^d x tq “ ^q l td 

Useful Load 

The useful loads applied to the alternator consist of series resistance and inductance 
as illustrated in figure 4. In terms of the direct and quadrature components of current, 
the load voltages are expressed as 

v d Bl Vid + X i<PW - X ^q 

v q* Vzq + X * (pl Zq ) + X J wi Zd 

These equations can be rewritten to eliminate the requirement of an explicit . 






Figure 4. - Schemati of useful load. 
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(This elimination simplifies the application of varied loads in the simulation. ) Dividing 
each side of the equation by gives 

v, 

— =hd + tan MpW " wtan h hq 

R i 


7T = i lq + tan + w tan 9 l Hd 

K l 

For a given power factor useful load, only need be varied to vary the magnitude of 
the load. 

The currents i^ and i ^ are components of the total current components i^ and 
*tq are determined from the alternator equations. Since the useful load current 

and the parasitic load current are assumed to comprise the total alternator current, 

*td = l l d + x pd 

X tq = + X pq 

The power dissipated in the useful load is equal to 

V l = v d i f d + v q i Z q 

Speed Controller 

The parasitic load speed controller is a control feature which distinguishes this 
electrical generating system from conventional turboalternator electrical generating 
systems. Figure 5 is a functional block diagram of the controller. The controller sim- 
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Figure 5. - Block diagram of parasitically loading speed controller. 
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ulated is a proportional control: output load current is proportional to the deviation of 
the input signal frequency from a nominal value Oq. 

When the controller is operating, the frequency deviation is detected and measured 
by the frequency discriminator. The output of the discriminator, amplified, determines 
the magnitude of the alternator current through the parasitic load resistor. 

In existing developmental parasitic load speed controllers (SNAP-2, SNAP-8, Bray- 
ton Power Systems) the control of the current through the load resistor is performed by 
phase control with saturable reactors or silicon controlled rectifiers (refs. 3 to 5, and 9). 
Such a phase-controlled current is illustrated in figure 6. Varying the firing angle 



(fig. 6) changes the effective value of this nonsinusoidal current. But in addition to 
changing the value of the current, varying the firing angle also has the effect of changing 
the phase angle between the applied voltage and the fundamental component of the non- 
sinusoidal current. The current fundamental lags the voltage as in an inductive circuit. 
It is this fundamental component of parasitic load current which is the power producing 
component. Thus, in effect, the variation of parasitic power is the result of a variable 
inductance. This effect is fully analyzed in reference 10. 
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The computer simulation of the speed controller is based upon this variable induc- 
tive effect of phase-controlled current. 

Frequency disc riminator. - The frequency discriminator is represented in the sim- 
ulation as a first degree lag network whose transfer function (ref. 11) is 

k d if 


where 


S2 - S2 n 

Af = f - f n = > 0 

U 2n 

(The zero limitation is imposed because a frequency less than fg implies a negative 
useful load, a nonrealistic quantity. See fig. 1.) The steady-state characteristic of 
the discriminator is illustrated in figure 7(a) where i^g is the output current of the 
discriminator at a frequency fg. Thus 

4 d = l b + 

Power amplifier. - The power amplifier has also been represented as a first degree 
lag network (ref. 11) whose transfer function is 


r A P+ 1 

The steady-state characteristic of this amplifier is illustrated in figure 7(b); thus 

*A = *A + *A0 

where i^g is the amplifier output current when i^ = 0, and is such a value that 

*A0 = -K A i D0 

The combined steady-state characteristic for the discriminator and amplifier is as shown 
in figure 7(c). 
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(a) Characteristic of frequency discriminator. 


j A 



(b) Characteristic of power amplifier. 



(c) Combined characteristic of discriminator and amplifier. 

Figure 7. - Form of steady-state characteristics of simulated speed 
controller components. 
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Figure 8. - Schematic of equivalent parasitic load. 


Parasitic load . - In order to simulate the effect of current phase control, the para- 
sitic load impedance is represented as a series fixed resistance and a variable induc- 
tance as in figure 8. 

Calculation of parasitic power . - The computer is programmed to balance the input 
and output powers of the alternator. The power provided by the turbine is presumed to 
be constant and is so set. The power is assumed to be completely transferred across 
the alternator air gap and dissipated in the field, the armature windings, and the useful 
and parasitic loads; that is, 


P ln = P f + 


P a + 


P 7 + P_ 


(For convenience P^ and P a may be considered part of P^ . ) 

In demonstrating the performance of the generating system a purposeful change is 
made in the useful load, producing a load torque change in the system. The resulting 
turbine and load torque unbalance produces a speed change. The rotary speed becomes 


(JL> = 

P J 

The parasitic load speed controller reacts to the speed change until the parasitic load 
change compensates for the useful load change and the turbine and load torques balance. 


12 



The power dissipated in the parasitic load resistor is 


P 


P 


= ij?R 

P P 


Since Rp is fixed, in order that the current indicated in this equation equal the para- 
sitic current generated by the speed control amplifier, the inductance value must vary. 
The analytical expression used to compute the reactance required to balance the system 
is 


X 


- v xp . ( v t - ‘p R p) 


1/2 


0)1 


O) 1 


In terms of the direct and quadrature components, the parasitic current becomes 


R 


R 


i =J-iv„->ii 


- u- i. 


~ pq x p p ' q x p p ~ pq p pd 


and the power dissipated in the parasitic load resistor is 


P p = Vpd + v q i pq 


Voltage Regulator-Exciter 


The model used for the voltage regulator is represented by a lead-lag network 
whose transfer function is given as 


v f = 


K(1 + Tl p) 
1 + TjP 


where the error signal 
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v e = v Ref - K l v t w 

This model was selected as a representative voltage regulator because it is the 
voltage control scheme used in an operating developmental generating system, the 
400-hertz Brayton Power System (ref. 12). A block diagram illustrating such a regu- 
lator is shown in figure 9. It represents a stabilized shunt field exciter, sensitive to 
changes in generator terminal voltage and frequency. 



Figure 9. - Simulated voltage regulator-exciter block diagram. 


COMPUTER REPRESENTATION OF GENERATING SYSTEM 
Computer Equations 

The equations presented in the preceding discussion were scaled and converted to 

a form convenient for application to a PACE 23 1R analog computer. The simulation in- 
volved 82 amplifiers, 32 multipliers, and two function generators. 

The parameter values used in the simulation are tabulated in table I. The values of 
this table were obtained from component test results - where available (refs. 13 and 14) - 
design goals (ref. 15), and estimates of the parameters of a 400-hertz developmental 
Brayton Power System. 

Block diagrams of the complete generating system as simulated are included as 
figures 10 to 12. 


Saturation Approximation 

The nonlinearity of the flux linkage function due to saturation can be approximated 
by a scheme discussed in reference 16. This scheme approximates the saturation by 
means of the assumptions that (1) only the direct axis mutual reactance varies with flux, 
and (2) the direct axis mutual flux is reduced as a result of a fictitious "saturation cur- 
rent" i This current is described in the following paragraphs. 

The alternator no load saturation curve (which is a measurable characteristic) 
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TABLE I. - PARAMETER VALUES 21 


[These values were used in simulation runs for 
which transient traces are shown in fig. 14.] 


Alternator 

Total inertia, J 

7. 28 

Armature resistance, R a 

0.0183 

Field resistance, R^ 

0. 00224 

Direct-axis damper resistance, R kd 

0. 1450 

Quadrature- axis damper resistance, R^ 

0. 1450 

Armature leakage reactance, X aL 

0. 1280 

Field leakage reactance, X fL 

0. 5770 

Direct-axis damper reactance, X kd 

0. 1158 

Quadrature-axis damper reactance, X^ 

0. 1680 

Direct-axis mutual reactance, X md 

1. 0000 

Quadrature -axis mutual reactance, X mq 

0.6230 

Speed control 

Discriminator current at fQ, i^Q 

0. 9550X10 3 

Amplifier gain, 

1. 8500X10 3 

Discriminator gain, 

0. 0464x10'^ 

Parasitic load resistance, R p 

0. 5556 

Amplifier time constant, r A , sec 

0. 0310 

Discriminator time constant, sec 

0. 0050 

Voltage regulator-exciter 

Gain, 

0. 0975 

Gain, K 

73.9 

Reference voltage, v Ref 

0. 0996 

Time constant, tj* sec 

0.095 

Time constant, sec 

1. 540 

Base values 

Field current, i f Q, A 

2.6 

Field voltage, v^q, V 

5760 

Frequency, f, Hz 

400 

Phase current, i tQ , A 

41. 6 

Power, P InQ , kW 

5 

Phase voltage, v tQ , V 

120 

Time base, a 

100 


a Per unit values unless otherwise indicated. 








Figure 12. - Block diagram of useful and parasitic loads. 


shows that the terminal voltage is a nonlinear function of field current (fig. 13(a)); that 
is, 


v t0 = ^mdO = CWV] (i f> 

where x mc j(if) is the nonlinear factor. The assumption that under load the same func- 
tional relation exists between current and reactance makes 

^md = [ X md^3 ^f " Hd + W 

If the saturation curve under load is plotted as shown in figure 13(b), then 








(a) Alternator no load characteristic. 



/ 

/ 

/ 


s 


(b) Saturation curve replotted to show i $ . 

Figure 13. - Use of saturation curves to determine saturation 
current. 


where the nonlinearity is included in the term This term is called the satura- 

tion current i g . Its value can be determined from the saturation curve as follows: 



td 


W " 


^md 

X mdO 


Then 


^md' X mdO ^f " l td + *kd " i s ) 
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Figure 14. - Simulation transient characteristics resulting from step changes in useful load. Useful load 
power factor = 0. 7 (lagging). 



In the simulation this saturation effect is included by using a nonlinear function 
generator to generate i g as a function of ^ mc j* 


Simulation Performance 

In order to demonstrate the performance of the simulation, recorder traces of 
several sample characteristics resulting from step changes in useful load are presented 
in figure 14. 

The characteristics display the expected theoretical performance of parasitically 
loaded electrical power generating system. Upon the step removal of the useful load 
the terminal voltage increases suddenly, returning to a steady-state value with a trans- 
ient typical of a lead-lag regulator characteristic. The slight increase in the useful 
load power immediately following the step change is the result of the voltage pulse from 
unloading of the alternator. The frequency increases exponentially to a steady- state 
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value, and displays the characteristic of an over-damped system. Similarly, the para- 
sitic load power increases in an exponential manner governed by the inductive nature of 
the parasitic load controller. The total power load (which in this case includes the arm- 
ature and field losses), after its transient dip, returns to its fixed steady-state value. 

When this load change is reversed, that is, upon the step application of essentially 
a per unit useful load, the displayed characteristics change in the opposite direction. 
There is one significant difference because (in this instance) the added useful load is in- 
ductive (0. 7 lagging power factor). The additional impedance to current change is evi- 
dent in the transient characteristics. 

CONCLUDING REAAARKS 

The simulation developed herein permits dynamic analysis of parasitically loaded 
electrical power generating systems. Analysis of such regulated power systems without 
resort to computers is prohibitive because of the complexity of the equations describing 
the systems. Used properly the simulation can be a valuable tool in the analysis of con- 
templated designs or the improvement of existing designs. With the computer, steady- 
state and transient stability can be observed as system conditions and parameters are 
changed. 

The simulation described can be altered and expanded to include additional opera- 
tional factors (such as variable turbine power) and more refined models of the compo- 
nents. With the addition of load-sharing circuits the simulation can be used for the study 
of paralleled alternators. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1970, 

120-27. 
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APPENDIX - SYMBOLS 


(Unless otherwise indicated, all quantities are expressed as per unit. 


f frequency, Hz (also functional nota- 

tion) 

i current 

i T current (variable component) 

J total inertia 

K gain 

P power 

p differential operator d/dt 

Q reactive volt-amperes 

R resistance 

T torque 

v voltage 

X reactance (inductive) 

a time base 

0 power factor angle 

r time constant 

rp flux linkage 

S2 electrical frequency, rad/sec 

w electrical frequency 


) 

Subscripts: 

A amplifier 

a armature 

D discriminator 

d direct axis 

e error 

f field 

In input 

k damper 

L leakage 

l useful load 

m mutual 

p parasitic load 

q quadrature axis 

Ref reference 

s saturation 

T turbine 

t total or terminal 

0 (zero) base or nominal value 
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